b=[19.97498333 39.94996666 19.97498333];
a=[68.2660915579 -0.0787823638 11.7126256859];
[h,w]=freqz(b,a,30)
plot(w/2/pi,abs(h)); grid on
disp(['w'])
(w/2/pi)'
abs(h)'